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ABSTRACT 



A rectangular cross section channel with both straight 
and curved test sections was used to investigate the effect 
of Tay lor-Goertler vortices on the heat transfer from the 
curved section by comparing the results with those obtained 
in the straight section. The flow was heated through the 
Joulean heating of one of the walls which was made of 
Temsheet. Liquid crystals were applied on the surfaces of 
both test sections. Experiments were run for several values 
of the volumetric flow rate corresponding to Reynolds num- 
bers between 258 and 920. The average Nusselt number was 
calculated for each experiment and plotted as function of 
the Reynolds and Dean numbers. Three distinct regions of 
the average Nusselt number were observed in the curved 
section when compared with the straight section. The presence 
of Taylor-Goertler vortices in the curved section was veri- 
fied from the distribution of the liquid crystal color 
bands. It was concluded that the overall effect of the 
vortices whenever they occurred was to improve the heat 
transfer process from the curved section. 
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I. 



INTRODUCTION 



A. THE TAYLOR-GOERTLER VORTICES 

Several studies have shown that a fully developed 
laminar flow of a viscous fluid along a concave wall does 
not remain two-dimensional. Instead, the flow forms a 
system of counter-rotating vortices - the Taylor-Goertler 
vortices - whose axes are in the direction of the main flow 
and with secondary velocities in both perpendicular direc- 
tions as shown in Figure 1. These vortices are caused by 
the variation in centrifugal forces on the fluid particles 
at different locations in the flow. 

In a curved channel such as used in the present inves- 
tigation the velocity profile of the fully developed laminar 
flow is nearly parabolic. The fluid near the center of the 
channel has a relatively higher velocity and is subjected 
to greater centrifugal forces than the fluid moving near 
the walls. Thus, the tendency is for the fluid in the cen- 
ter of the channel to move outward toward the concave wall. 
The fluid near the concave wall is unable to resist this 
action and must move in the spanwise direction and then 
radially inward replacing the central fluid. Once in the 
center of the channel, this fluid is subjected to the higher 
flow velocity and the rotation is continued. This cyclic 
motion forms the counter-rotating Taylor-Goertler vortices. 
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Figure 1. Sketch of Taylor-Goertler 
Vortices 
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Although the Taylor-Goertler vortices are laminar 
vortices, they may affect the transition from laminar to 
turbulent flow and are therefore important to understanding 
this phenomenon. 

The influence of the Taylor-Goertler vortices may also 
explain the increase in the observed heat transfer rates 
from curved walls. 

The Taylor-Goertler vortices have similarities to other 
vortex flows such as the longitudinal vortex rolls developed 
in the laminar forced convection heating of fluid layers 
between parallel plates. The striations seen at stagnation 
points on bluff bodies and the cross hatching observed in 
reentry vehicles have also been explained, at least in 
part, by the presence of Taylor-Goertler vortices. 

There are many possible applications that could result 
from investigation and understanding of the Taylor-Goertler 
vortices. A few of such applications include improved 
cooling of turbine blades and other external surfaces, 
reduction of pressure losses in bends and a better under- 
standing of the transition from laminar to turbulent flow. 
Heat exchangers could be designed to take advantage of the 
low pressure drop associated with laminar flow and the 
improved heat transfer characteristics of the flow in the 
presence of the vortices. 
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B. BRIEF HISTORY 



Lord Rayleigh [Ref. 1] in 1916 was the first to consider 
the instability of an inviscid fluid rotating symmetrically 
about an axis. In his study he determined the stability 
criterion for the inviscid fluid to be that the product 
of the local circumferential velocity and' the local radius 
of curvature either increases or at least remains constant 
as the radius increases. 

In 1923 G.I. Taylor [Ref. 2] in an extensive analytic 
and experimental study extended the work of Rayleigh to 
viscous fluids. In his investigation for the flow between 
two rotating cylinders, he found that when the inner cylin- 
der rotates and the outer cylinder is stationary the motion 
of the fluid becomes unstable when the value of the Taylor 
number defined as: 



exceeds 41. In the equation above, d is the spacing between 
the cylinders, assumed small when compared with Ri , the 
radius of inner cylinder and Re is the Reynolds number based 
on the circumferential velocity of the inner cylinder and 
the length d. For values of the Taylor number greater than 
the critical value a secondary motion sets in and the 
Taylor vortices occur. 
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A similar type of instability is observed when a 
viscous fluid flows in a curved channel due to a pressure 
gradient acting around the channel. This problem was first 
considered analytically by W.R. Dean [Ref. 3] for a channel 
formed by two concentric cylinders where the spacing between 
the cylinders was small when compared with the radius of 
the inner cylinder. Dean concluded that the instability 
will first arise and secondary flow vortices similar to the 
Taylor vortices will form when the Dean number defined as 



is greater than 36. In this case the Reynolds number was 
based on the mean velocity of the unperturbed flow and the 
length d. 

The analytical work of Dean was later verified by W.H. 

Reid [Ref. 4] using a simplified method. H. Goertler [Ref. 

5] in 1940 studied the stability of laminar boundary layer 
profiles on slightly curved walls relative to small distur- 
bances. Goertler found that these disturbances were similar 
to those investigated by G.I. Taylor which led to the 
instability in the form of vortices. From numerical calcu- 
lations he concluded that amplified disturbances were pro- 
duced only on the concave walls. The approximate results 
obtained by Goertler were verified in 1955 with an exact 
solution developed by G. Hammerlin as reported by H. Schlichting 
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[Ref. 6] and more recently by A.M.O. Smith [Ref. 7] through 
an extensive numerical analysis. 

In 1955, F. Kreith [Ref. 8] investigated the influence 
of curvature on heat transfer for fully developed turbulent 
flows and reported that the heat transfer from a heated 
concave wall was considerably higher than from a convex 
wall of the same curvature and under the same conditions 
of flow. 

L. Persen [Ref. 9] in 1965 was able to relate the increase 
in the heat transfer rate from curved walls with the presence 
of the Taylor-Goertler vortices. Persen considered the 
special cases of a very high and very low Prandtl number and 
determined that the overall effect of the vortices was to 
increase the rate of heat transfer through the boundary 
layer. 

In 1970 P. McCormack, et al. [Ref. 10] reported the first 
experimental work on the effect of the Taylor-Goertler vor- 
tices on the heat transfer. It was pointed out in their 
study that in order to explain theoretically the effect of 
the vortices, it was necessary to retain the non-linear 
terms in the flow equations. This was considered by R. 
Kahawita and R. Meroney [Ref. 11] who concluded that the 
high-order terms and the normal velocity components of the 
main flow become increasingly important at small wave 
numbers . 
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Y. Mori and Y. Uchida [Ref. 12] in 1966, considered 
the fully developed forced convective heat transfer between 
horizontal flat plates. They observed that when the tem- 
perature difference between the plates increased above a 
critical value, there was a formation of longitudinal vortex 
rolls, similar to the Taylor-Goertler vortices, with the 
axes parallel to the flow direction. M. Akiyama, et al., 
[Ref. 13] in 1971 confirmed the onset of these vortices at 
a critical Rayleigh number of 1708 for the case of heating 
from below. 

The published literature concerning the flow and heat 
transfer in curved rectangular channels has been rather 
limited as reported by R. Shah and A. London [Ref. 14]. 

In 1970, K. Cheng and M. Akiyama [Ref. 15] developed a 
numerical solution for the forced convection heat transfer 
in a curved rectangular channel but only for small aspect 
ratios (nearly square cross section) . Analytic and experi- 
mental results for a fully developed constant wall heat 
flux square cross section channel were obtained by Y. Mori, 
et al. [Ref. 16] . K. Cheng, et al. , [Ref. 17] in 1974 
solved the Graetz problem in a curved square channel using 
a numerical method. 
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II. INTENT OF THE STUDY 



The purpose of this study was to investigate the 
effects of the Taylor-Goertler vortices on the heat trans- 
fer in a curved rectangular channel and to compare with the 
heat transfer in a straight channel of similar cross section. 
It was expected that the effect of the vortices would result 
from the secondary velocity components transporting heated 
fluid from the curved wall toward the center of the channel 
and carrying cooler fluid from the center to the curved 
wall. It was also expected that the heat transfer process 
from the curved wall would show improvement as a consequence 
of the presence of the vortices. 

For the purpose of this study, a rectangular cross 
section channel with both straight and curved test sections 
was considered. The results obtained in the curved section 
were compared with those obtained in the straight section 
to determine the effect of the Taylor-Goertler vortices on 
the heat transfer. 

An analytical solution for the case of the straight 
test section was given by P. McCuen, et al., [Ref. 18] in 
their 1962 study. The results obtained in the straight 
section were compared with the analytical results of McCuen. 

For the curved section, the results obtained in the 
present investigation were compared with the numerical 
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solution developed by K. Cheng and M. Akiyama [Ref. 15] 
for the heat transfer in rectangular curved channels with 
small aspect ratios. 
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III. EXPERIMENTAL WORK 



A. DESCRIPTION OF THE APPARATUS 

A rectangular cross section channel built as described 
in Refs. 19 and 20 was modified to include two heated 
sections in order to achieve the objective of this study. 

The channel was made of two 0.635 cm thick sheets of Plexi- 
glas separated by 0.635 cm thick spacers which also formed 
the sides of the channel. 

The channel was composed of a straight portion 122.0 
cm long, a curved section which completed a 180 degrees 
turn and a final short straight section as shown in Figure 
2. The radius of curvature of the interior concave wall 
of the curved portion was 30.5 cm. The cross section of the 
channel was 0.635 cm high and 25.4 cm wide with an aspect 
ratio of 40 and a cross sectional area of 16.129 square 
centimeters. The cross section of the channel at both test 
sections is shown in Figure 3. 

The working fluid was air which entered the channel 
through an aluminum entrance nozzle attached to the 
beginning of the straight portion and which was covered 
with cheesecloth during the experimental work. At the end 
of the channel, the flow was conducted through an aluminum 
exhaust nozzle connected by flexible tubing to a Fisher and 
Porter Company variable area flow meter, model number 
10A3565. This rotometer had a 100% full scale flow rate 
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Figure 2. Schematic Illustration of Apparatus 
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Figure 3. Cross Section of Channel at Both Test Sections 





of 0.314 cubic meters of air per minute (11.1 standard 
cubic feet of air per minute) . 

The flow of air was drawn through the channel, the 
flexible tubing and the rotometer by two electrically driven 
Cadillac centrifugal blowers, model G-12. The blowers 
were connected in parallel, and their speeds were controlled 
by varying the motor voltage with a General Radio Company 
Variac Autotransformer type W10MT3. The voltage was 
regulated by a Sorenson A.C. voltage regulator, model 
R1050. 

To obtain the experimental heat transfer data, two 
test sections were constructed: a straight test section 

29.2 cm long in the straight portion of the channel after 
the hydrodynamic entrance region and a curved test section 

28.3 cm long subtending an arc of 53.1° in the curved 
portion of the channel. The area of the straight test 
section was 741.93 square centimeters and that of the curved 
test section was 717.74 square centimeters. 

In each test section the outer wall was removed and 
modified as shown in Figures 4 and 5. A piece of Temsheet - 
a carbon impregnated porous paper with the property of 
uniform electrical resistivity - was glued to the entire 
interior surface of the outer wall in each section. The 
flow of air was heated through the Joulean heating of the 
Temsheet. As the electrical resistance of the Temsheet was 
not constant but slowly variable with the temperature, a 
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Plexiglas 



Thermocouple 
suspended 
in the flow 




Figure 5. Internal Arrangement of the 
Curved Test Section 
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precision resistor with an electrical resistance of 2.013 
ohms was inserted into the heating circuit in series with 
the Temsheet to calculate the instantaneous power supplied. 

The electrical current to heat the Temsheet was supplied 
from a LAMBDA Regulated Power Supply model LK345A FM. 

The variables to be measured were: the flow inlet 

temperature at the entrance of the channel (T ) / the flow 

outlet temperature after each test section ( T 0UT ) / the wall 

temperature of the heated plate in each test section ( T^ ) , 

the flow exit temperature at the end of the channel (T t?vt _) , 

ilXJ- 1 

the temperature of the insulation for each test section 
(Tj^s) an< 3 the voltages across the precision resistor (V p ^) 
and the heater (V ) . 

n 

Thermocouples were used to measure the temperatures. 
Cholesteric liquid crystals were used also to measure quali- 
tatively the surface temperature distribution on the Temsheet. 
Liquid crystals exhibit dramatic changes in color for small 
changes in temperature. (An excellent reference on the use 
of liquid crystals is the article by Cooper, et al. [Ref. 

21]). The liquid crystals would provide a reasonable 
qualitative comparison between the heat transfer process 
on both test sections through the observation of the dis- 
tribution of the color bands in each case. The occurrence 
of color stripes on the curved wall, corresponding to areas 
of relatively high and low temperature as a result of the 
presence of the Taylor-Goertler vortices was expected to 
be observed. 
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To select the appropriate liquid crystals to be applied 
on the surface of the Temsheet, the values of the flow inlet 
temperature, the Reynolds number, the Prandtl number and 
the Nusselt number were assumed and an approximate energy 
balance was established: 

Q = m Cp (t qut - t in ) = F A pL at 



or 



AT 

T - T 
OUT IN 



Re Pr , Ac . 
Nu A pL 



where 



AT 



Two 




+ T 

OUT 

2 



For a desired difference between the outlet and the 
inlet temperatures of fifteen degrees Celsius, a wall 
temperature of fifty-five degrees Celsius was calculated. 

The required power considering the possible losses was 
estimated between 45 and 75 watts. 

Based on the calculated wall temperature the NCR Micro- 
encapsulated Liquid Crystal R-53 (onset of red at fifty 
three degrees Celsius) was selected. After some preliminary 
experiments, it was found that the response curve of the 
liquid crystal R-53 allowed for a large variation of the 
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wall temperature within the range of the color blue. It 
was then decided to use a combination of three different 
liquid crystals instead of a single one in order to obtain 
several color bands and a better appreciation of the average 
wall temperature. A new Temsheet test section was constructed 
and covered with a mixture of the NCR Microencapsulated 
Liquid Crystals R-49, R-56 and S-62 prepared in accordance 
with Ref. 21 and applied directly to the surface with a 
spray gun. About twenty coats were applied to each plate 
to obtain a clearly observable differentiation between the 
different color bands which it was hoped could be recorded 
photographically . 

Thirty seven copper-constantan 30 gauge thermocouples 
were constructed and inserted at several locations in the 
apparatus to measure the different required temperatures. 

All the thermocouples were connected to a Honeywell thermo- 
couple switch mounted on the side of the apparatus. (See 
Figure 6 . ) In each test section five thermocouples were 
inserted between the first and the second layers of the 
insulation and another five between the second and third 
layers. These thermocouples were connected in parallel 
to read an average value of the temperature of the insula- 
tion (T ) . Twelve thermocouples, individually connected, 
were suspended in the flow, four in the entrance of the 
channel, four at the end of the straight test section and 
the last four at the end of the curved test section. At 
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each one of these locations the thermocouples were distributed 
spanwise at intervals of 5.0 cm. These thermocouples were 
used to measure the bulk temperature of the air. In each 
test section two thermocouples were placed in direct con- 
tact with the Temsheet through holes of very small diameter 
drilled on the middle of each plate at a distance of 10.2 
cm from the leading and trailing edges. The beads of these 
thermocouples were electrically insulated with ENMAR Heat 
Resisting Glyceryl Phthalate. A final thermocouple was 
inserted into the flexible tubing close to the rotometer. 

All the thermocouples were calibrated using a ROSEMOUNT 
Commutating Bridge, model 920A, and a ROSEMOUNT Calibration 
Bath, model 913A. 

The thermal insulation consisted of ARMSTRONG ARMAFLEX 
22 Sheet Insulation - a flexible foamed plastic material. 

Three layers of insulation, each one with a thickness of 
0.635 cm and covering an area slightly larger than each 
test section, were used. The insulation was fixed in place 
with the use of Adhesive Corrosion Resistant tape as shown 
in Figure 7. 

To prevent any leakage of air into the channel, the 
whole apparatus was sealed with DOW CORNING 781 building 
sealant. 

Two aluminum electrodes 0.318 cm thick were inserted 
between the Temsheet and the plexiglas in each test section 
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Figure 7. General Arrangement of Apparatus and Instrumentation 
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as shown in Figure 3. The plates with the Temsheet, the 
liquid crystals and the aluminum electrodes were assembled 
in the channel with plastic screws for the purpose of 
electrical insulation. The Thermocouple Switch was con- 
nected to a DYMEC Integrating Digital Voltemeter model 
2401B, and to a reference ice bath. A HEWLETT PACKARD 
Digital Recorder, model 562A, provided an automatic printed 
record of the instantaneous data for each position of the 
thermocouple switch- (See Figure 8.) 

The general arrangement of the apparatus and the 
instrumentation is sketched in Figure 2. 

B. EXPERIMENTAL PROCEDURES 

Experiments were run in each test section for several 
values of the volumetric flow rate in the range between 
0.063 and 0.220 cubic meters of air per minute. 

The appropriate electrical power required for each 
value of the volumetric flow rate in order to heat the 
Temsheet up to a temperature at which the liquid crystals 
would change colors was not known. Several preliminary 
experiments were performed to determine those values of 
electrical power and also to obtain an appreciation for 
the general development of the experiments. 

Estimates of the required time to reach steady-state, 
the necessary data to be taken, the frequency of the 
readings, the overall time of each experiment and the 
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interval between experiments were obtained from the 
preliminary experiments. 

The experiment was considered to have reached steady- 
state when the variation in the average wall temperature 
given by the thermocouples in contact with the Temsheet 
was less than one degree Celsius for a period of ten 
minutes . 

The time to attain steady-state was found to be between 

one hour and one and one half hours. It was decided that 

the duration of all experiments would be two hours. Data 

* 

was taken each half hour during the first hour and every 
ten minutes in the second hour. 

As the area of the curved test section was slightly 
less than the area of the straight test section, the instan 
taneous power supplied was different. In order to Obtain 
similar values of dissipated power in both test sections 
for each value of the volumetric flow rate, the power had 
to be adjusted for the curved section during the experiment 

The room temperature was determined with a mercury 
thermometer located in the vicinity of the apparatus. The 
volumetric flow rate was read directly from the rotometer. 
The instantaneous values of Tw, T EXIT , T INg , T JN and T QUT 
were automatically printed by the digital recorder in 
millivolts. Based on the results of the calibration of 
the thermocouples, the standard thermocouple tables were 
used to convert millivolts to degrees Celsius. 
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The voltages across the precision resistor (V__) and 

P R 

the heater (V^) were recorded in volts, and the instantaneous 
power supplied (Qp) was then computed using the expression: 



V 



PR 



V 

r pr h 



where Rpp was the electrical resistance of the precision 
resistor, previously determined as 2.013 ohms. 

Prior to each experiment, an ice reference bath was 
prepared and used for the single reference thermocouple 
connected in series with the thermocouple switch. 

The steady-state distribution of the liquid crystal 
color bands on the surface of both test sections was 
sketched for most of the experiments and photographed for 
the case of the curved test section using a Nikomat 35 mm 
camera with a Nikon 55 mm f 3.5 Macro-Nikkor lens. 

During the preliminary experiments the bulk temperature 
thermocouples produced anomolous readings. Determining 
the cause required considerable time. An energy balance 
on the straight test section indicated the value of power 
transmitted by the heater to the air flow to be greater 
than the power supplied to the heater. 

It was thought at the beginning that leakage of air 
into the channel could be the cause. However it was later 
shown that the problem consisted of erroneous temperature 
readings given by the thermocouples suspended in the flow. 
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These thermocouples had been inserted into the channel 
through Plexiglas plugs which could be removed for replace- 
ment of any defective thermocouple. 

Several tests were performed, and it was determined 
that there was an electrostatic interaction between the 
Plexiglas plugs and the beads of the thermocouples. The 
recorded values of the temperature were from two to five 
degrees Celsius higher than the actual values. This dis- 
tortion was more evident for the higher values of the 
volumetric flow rate. The response time of the thermo- 
couples also was much slower than normal. This was con- 
firmed by removing the thermocouples from the apparatus and 
observing the approach to room temperature. It was found, 
however, that once the thermocouple beads were removed 
from the vicinity of the Plexiglas plugs they would read 
correct values of temperature with a normal response time. 
As the height of the channel was only 0.635 cm, it was 
difficult to move the thermocouple beads away from the 
Plexiglas plugs. 

Considering that the problem probably resulted from 
some electrostatic charge induced in the Plexiglas and 
transmitted to the beads it was decided to try different 
materials for the construction of the plugs. Phenolic 
resin fiberboard and Dine wood were used, but the same 
unsatisfactory results were obtained. 
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Finally it was decided to devise a process to insert 
the thermocouples into the channel without the use of the 
plugs and also to bend the thermocouples in order to place 
them parallel to the air flow as shown in Figures 4 and 5. 
When this was accomplished, the correct data was obtained. 
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IV. PRESENTATION OF DATA 



A. ANALYSIS 

Since Temsheet has the property of uniform electrical 
resistivity, the heated wall was considered to be a constant 
heat flux surface. Also because of the large aspect ratio, 
the channel was considered to be approximately parallel plates. 
It was therefore assumed that the experimental configuration 
represented forced convection between parallel plates with 
one surface subjected to constant heat flux and with some 
heat loss at the opposite surface. To analyse this situation 
several quantities were defined and calculated as follows: 

The average heat transfer coefficient between the heated 
wall and the flow of air in each test section was defined 
by the equation 

Qair = h A pL AT 

In this equation, Qair was the heat convected- (See the 

Figure in Appendix B.) h was the average heat transfer 

coefficient, A was the area of the heated surface and 
P Li 

AT was the difference between the average wall temperature 
(Two) and the bulk temperature of the fluid (T p ) . The bulk 
temperature was defined as the arithmetic average of the 
fluid inlet temperature (T IN ) and the fluid outlet temperature 
( T OUT ) ’ In mathematical form AT was expressed by: 
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T 



IN 



+ T 



OUT 



AT 



Two - T 



B 



Two 



2 



The heat convected (Qair) was calculated with the 
equation 



where Cp was the specific heat of air at constant pressure 
and m was the mass flow rate of air given by the product 
of the volumetric flow rate (q) and the density of the air 
(p) . This density was computed assuming perfect gas 
behavior at atmospheric pressure and using the flow exit 
temperature (T EXIT ) • 

The average Nusselt number was then computed from the 
average heat transfer coefficient through the expression: 



where d was the height of the channel and Kair was the 
thermal conductivity of the air at the bulk temperature 



Qair 





The Reynolds number was calculated for all the 
experiments and defined as: 



pUd 
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where 



U = 



_q_ 

Ac 



Ac was the cross sectional area of the channel and y was 
the dynamic viscosity of the air at the flow exit temperature 

(T EXIT ) ' 

For the experiments in the curved test section, the 
Dean number defined by: 

De = Re 




was also calculated. In this equation, d was the height 
of the channel and Ri was the radius of curvature of the 
interior surface of the convex wall. 

The heat lost by conduction through the outer plate 
and the insulation was determined with the expression: 



Lo 



iT INS 

(iX INS /K INS A PL> 



where AX . TC was the thickness of each layer of insulation, 
K INS was tiie therma l conductivity of the insulation and 
AT was the difference in the temperatures given by the 

_L IN O 

thermocouples located between the first and second layers 
and those located between the second and third layers. 
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In order to verify that Qair as measured by the 
temperature rise in the air was actually the convection 
heat transfer from the Temsheet, an estimate of the radiation 
heat transfer in the channel had to be made. To facilitate 
this, one experiment was run in which the temperature of 
the lower plate was measured with two thermocouples. For 
this experiment the radiated heat transfer was computed 
with the expression: 



4 4 

n _ a (Two - Twi ) 

R - r r 

where a is the Stef an-Boltzmann constant and R was the 

K 

total resistance to radiation heat transfer between 
the plates given by: 



R 



R 





1 ] 



In this expression is the emissivity of the interior 

surface of the outer plate (Temsheet) and £ ^ is the 
emissivity of the interior surface of the inner plate 
(Plexiglas ) . 

A sketch of the control volume for the energy balance 
indicating the different components involved and the sample 
calculations, for the case in which the radiated heat 
transfer was also computed, are given in Appendix B. 



43 



u 

tH AV 






B. RESULTS 



The data collected from the experiments was evaluated 
through the expressions described in the previous section 
as illustrated in Appendix B, and the results for the major 
variables involved are shown in Table I for the straight 
test section and in Table II for the curved test section. 

A plot of the natural logarithm of the Reynolds number 
versus the natural logarithm of the average Nusselt number 
is given in Figure 9 for both test sections. The incertainty 
bands are shown for two representative points. Sample 
calculations for obtaining these incertainty bands are 
shown in Appendix A. Figure 10 is a plot of the natural 
logarithm of the Dean number versus the natural logarithm 
of the average Nusselt number for the curved test section. 

The results indicated an increase in the heat transfer 
rates with increasing values of Reynolds number for both 
test sections. However, when comparing the plotted values 
of the average Nusselt number for the curved section with 
those of the straight section, three distinct regions were 
detected in the curved section. For the highest values 
of Reynolds number corresponding to the highest values of 
the volumetric flow rate, the Nusselt numbers for the 
curved test section were greater than those of the straight 
test section. For the lowest values of Reynolds number, 
the opposite was noticed with greater heat transfer rates 
occurring in the straight test section. A transition region 
was observed for intermediate values of Reynolds number. 
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Summary of Results for Curved Section 
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Figure 9. Natural Logarithm of Reynolds Number 
versus Natural Logarithm of Average 
Nusselt Number for Both Test Sections 
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Figure 10. Natural Logarithm of Dean Number versus 
Natural Logarithm of Average Nusselt 
Number for Curved Test Section 
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These differences in heat transfer behavior for the 
curved section as compared with the straight section were 
then related to the occurrence of Taylor-Goertler vortices 
through the observation of the distribution of the liquid 
crystal color bands on the surface of the curved section 
for the different values of Reynolds number. It was observed 
that for the smallest values of Reynolds number the differ- 
ent color bands occurred in almost straight lines perpendicu- 
lar to the direction of the flow. The same had been verified 
in the straight section independently of the value of Reynolds 
number. It was then observed that as the Reynolds number 
was increased for the curved section, color stripes oriented 
in the streamwise direction, became visible in the region 
close to the trailing edge. These color stripes were ex- 
plained by the presence of the Taylor-Goertler vortices 
and, for the highest values of the Reynolds number, could 
be seen over the entire surface of the curved test section. 
Schematic drawings of the distribution of the color bands 
on the surface of the straight and curved test sections are 
shown in Figures 11, 12, 13 and 14. 

The occurrence of color stripes and consequently Taylor- 
Goertler vortices at the same values of the Reynolds number 
for which the average Nusselt numbers were greater in the 
curved test section led to the conclusion that the effect 
of the vortices was to increase the heat transfer rates 
from the curved section. 
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Figure 11. Distribution of Liquid Crystal Color Bands 
in Straight Test Section for Re = 526 
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Figure 12. Distribution of Liquid Crystal Color Bands 
in Curved Test Section for Re = 638 (De = 
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Figure 13. Distribution of Liquid Crystal Color Bands 
in Curved Test Section for Re = 386 (De = 



♦ 
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Figure 14. Distribution of Liquid Crystal Color Bands 
in Curved Test Section for Re = 258 (De = 



The plot of the natural logarithm of the Dean number 
versus the natural logarithm of the Nusselt number for the 
curved section was then analyzed, and straight lines were 
fitted through the results in each region using least 
squares regression. These are shown in Figure 10. The 
slopes of these lines for the regions of highest and lowest 
values of Dean number were determined as respectively 0.557 
and 0.591 while in the transition region the slope was 
1.57. The limiting values of Dean number for the transition 
region were calculated from the plot as 57.1 and 67.4. 

Three correlations expressing the Nusselt number as a func- 
tion of the Dean number for each one of the regions were 
determined as follows: 

1. Low Values of Dean Number (De < 57.1) 

Nu = 0.292 De 0 * 591 

2. Transition Region (57.1 < De < 67.4) 

Nu = 0.00556 De 1 * 57 

3. High Values of Dean Number (De > 67.4) 

Nu = 0.395 De 0,557 
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A single correlation for the overall range of Dean 
numbers was also determined as: 



Nu 



0.171 De° 



.739 
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V. DISCUSSION AND CONCLUSIONS 



The measurement of the surface temperature of the 
unheated plate in the straight test section showed that 
this temperature was very close to the bulk temperature 
of the fluid. 

It was then assumed that in the absence of significant 
difference in temperature between the inner wall and the air 
flow, the heat convected to the air was exclusively from 
the heated wall. Based on this assumption it was decided 
to compare the results of the present study for the straight 
test section with the analytical solution of the same prob- 
lem developed by McCuen, et al. [Ref. 18] . The hydraulic 
diameter of the channel was calculated as Dh = 1.239 cm, 
and the values of Reynolds number based on the hydraulic 
diameter as used by McCuen also were determined. A dimen- 
sionless axial length coordinate x, given by 

— _ x 

x “ Dh Re Pr 

where Pr was the Prandtl number, was calculated for 
x = 29.2 cm (length of the straight test section) for all 
the experiments. Based on the values of x it was determined 
that the development of the experiments would occur essen- 
tially in the thermal entrance region, and it was decided 
to integrate the expression of the local Nusselt number 
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for small values of x 



Nu 



0.67095 x 1/3 



_JL 

- 2(1 



Qwi, 

Qwo' 



X 



given by McCuen, et al. , in order to obtain the following 
expression for the average Nusselt number. 



Nu = - 



4(1 + 



Qwo ; 



In [1 - 



0.67095 



(1 + 



Qwi,- 2 / 3 , 

Qwo ; J 



In this expression Qwi was the heat convected to the air 
from the unheated wall, and Qwo was the heat convected 
from the outer Temsheet wall. The average Nusselt numbers 
were then calculated with the expression developed above 
where, in accordance with the previous assumption, the 
value of Qwi was taken as zero and the value of Qwo was 
considered to be the same as Qair. 

The average Nusselt numbers obtained from the experi- 
mental data in the straight section were recalculated as 
function of the hydraulic diameter and were plotted versus 
Reynolds number as shown in Figure 15 along with the 
analytical expression from McCuen. The only significant 
difference observed for the lowest value of Reynolds 
number may be explained by the fact that the value of x 
was equal to 0.064. In this case the Nusselt number 
expression for small values of x is probably not valid. 
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Figure 15. Comparison of Present Data with the Study by McCuen, et al. 



A comparison of the data from the curved section with 
the analytical study by Cheng and Akiyama [Ref. 15] was 
attempted also for the curved test section, but the results 
were less conclusive. Their investigation was done for 
small aspect ratios (nearly square cross section) , and the 
results were plotted only for small values of Dean number. 
An extrapolation of these results to greater values of Dean 
number for the case of the largest aspect ratio (5 to 1) is 
shown in Figure 16 simultaneously with the results of the 
present study. 

After these comparisons, several conclusions about the 
heat transfer in rectangular channels were developed from 
the results obtained in the present investigation. The 
heat transfer rates were demonstrated to increase with 
increasing values of Reynolds number in both test sections. 
The values of the average Nusselt number for the curved 
section when compared with the values for the straight 
section seemed to indicate that there are three distinct 
regions of Reynolds or Dean number which can be considered. 
These regions were: an initial region for the lowest 

values of Dean number where the heat transfer rates from 
the straight test section were greater than from the 
curved test section, a final region for the highest values 
of Dean number where the opposite was true and a transition 
region for the intermediate values of Dean number. 
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Figure 16. Comparison of Present Data with the Study by 
K. Cheng and M. Akiyama 



The observation of the presence of Taylor-Goertler 
vortices at values of Dean number for which the heat 
transfer rates were greater from the curved test section 
led to the conclusion that the overall effect of the 
vortices was to improve the heat transfer process from the 
curved section whenever their presence occurred. The 
critical value of Dean number for the onset of the vortices 
was determined from Figure 10 as 57.1. 

Finally it seemed apparent from the comparison with 
the work of Cheng and Akiyama that for values of aspect 
ratio greater than five the results are independent of 
the aspect ratio. 
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VI. RECOMMENDATIONS 



Several recommendations can be made to improve the 
experimental study of heat transfer in curved as well as 
straight rectangular channels. More experiments should 
be performed in the transition region of the Dean number 
for the curved test section. This would permit a better 
understanding of the variation of the average Nusselt number 
with the Dean number and a more rigorous estimate of the 
critical value of Dean number above which the Taylor- 
Goertler vortices will form. The effect of the vortices 
on the heat transfer from the curved section could then be 
explained with more detail, and improved correlations for 
the average Nusselt number as a function of the Dean number 
could be developed. 

With respect to the apparatus, it became apparent that 

an improved method of suspending the thermocouples in the 

flow to measure T t „ and T^ rTm should be devised. Considerable 

IN OUT 

thought must be given to this problem in order to eliminate 
the interaction between the Plexiglas and the beads of the 
thermocouples and also to make the distance between the 
beads and the plates more uniform for all the suspended 
thermocouples. The surface temperature of the inner wall 
should be measured in future experiments. 

Finally, a motion picture study of the onset of the 
Taylor-Goertler vortices in the curved test section, using 
the liquid crystals, would be possible and of great interest. 
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APPENDIX A 



ERROR ANALYSIS 



The uncertainties for the major variables in the 
experiments were calculated in accordance with the method 
described by S. Kline and F. McClintoch [Ref. 22] . The 
estimates of the uncertainty in the measured quantities 
were made quite conservatively so that there was considerable 
confidence in the calculated uncertainties. As an example, 
the calculation of the uncertainty in the Reynolds number 
is given below. The Reynolds number was defined by the 
equation : 



of the uncertainties in the cross sectional area of the 
channel and in the reading of the volumetric flow rate. 

The uncertainties in the density, the height of the channel 
and the dynamic viscosity were obtained from estimates as 



Re 



pUd 

V 



and the uncertainty was calculated as: 



dRe 

Re 




The uncertainty in the flow velocity was determined as 
tr. = 0.02176 from similar calculations based on estimates 
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— = 0.00082, ^ = 0.02000 and as 0.00223. The 
P d y 

uncertainty in the Reynolds number was then calculated 
as 



dRe 

Re 



(0. 00082) z + (0. 02176) + (0.02)^+ (0.00223) 



0.02965 



The uncertainty in the Nusselt number was computed 
with the expression: 



dNu 

Nu 



J f) > 2 + <^> 2 + <f ) 2 + o 2 + <^> 2 + # > 2 + < 



dT -T 
OUT 

T -T 
CUT 



IN) 2 
IN 



The values of the uncertainties for other variables 
are listed below. 



Quantity Uncertainty 



Ac 


0.02010 


A 


0.00265 


PL 




CP 


0.00415 


De 


0.03576 


h 


0.06242 


Kair 


0.00038 


Nu 


0.06555 


q 


0.00835 


Qair 


0.04050 


Re 


0.02965 
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Quantity 



Uncertainty 



D 


0.01640 


t in 


0.00453 


T 

OUT 


0.01060 


T 

W 


0.03207 


T -T 

OUT IN 


0.03940 


AT = T r7 - T_ 
W B 


0.04743 


y 


0.00223 


p 


0.00082 
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APPENDIX B 



SAMPLE CALCULATIONS 

A sketch of the control volume for the energy balance 
on the straight test section indicating the major heat 
transfer components involved is shown in Figure 17. A 
similar diagram could be considered for the curved test 
section. Sample calculations of the power supplied (Op) , 
the heat convected to the air (Oair) , the heat lost through 
the insulation (Q ) , the radiated heat transfer (Q ) , 

JjO k 

the average heat transfer coefficient (h) , the average 
Nusselt number (Nu) and the Reynolds number (Re) are given 
below for the straight test section. 
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ENERGY BALANCE 




Figure 17. Energy Balance in Straight Test Section 

SAMPLE CALCULATIONS 
A . DATA 



t room 


— 


23.8 0 


C 






q 


= 


0.126 


m'Vmin 






V PR 




3.223 


V 






V H 


= 


35.580 V 






r pr 


= 


2.013 


SI 






T 

INS1 


= 


1.997 


mv = 


49.13 


°C 


T 

INS2 


= 


1.542 


mv = 


38.33 


°C 


T 

INI 


= 


0.963 


mv = 


24.29 


°C 


T 

IN2 


= 


0.962 


mv = 


24.27 


°C 



I 
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T 

IN3 




0.955 mv = 


24.10 


°C 


T 

IN4 


= 


0.965 mv = 


24.34 


°C 


T 

EXIT 


= 


1.072 mv = 


26.97 


°C 


Twol 


= 


2.684 mv = 


64.93 


°C 


Two 2 


= 


2.432 mv = 


59.19 


°C 


t outi 


= 


1.591 mv = 


39.52 


°C 


T OUT2 


— 


1.606 mv = 


39.87 


°C 


T OUT3 


= 


1.527 mv = 


37.97 


°C 


T 0UT4 


= 


1.662 mv = 


41.21 


°C 


Twil 


= 


1.417 mv = 


35.34 


°C 


Twi2 


= 


1.265 mv = 


31.66 


°C 


a pl 


= 


0.0742 m 2 






A 

c 


=: 


0.0016 m 2 






d 


= 


0 . 00635 m 






Kair 


= 


0.02662 W/m 


°C 




Cpair 




1.0057 KJ/Kg 


°C (0 


.241 BTU/lbm 


pair 


= 


1.983 x 10 -5 


Kg/m* 


sec 


k ins 


= 


4.18 x 10" 2 


W/m °C 




ax ins 


= 


0.00635 m 






e 

wo 


=r 


0.90 






£ . 
W1 




0.75 






Fwo-wi = 


1.0 






0 


— 


5.669 x 10 -8 


W/m 2 ° 


K 
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B. TEMPERATURE CALCULATIONS 

1- Average Inlet Temperature (T ) 

T + T 4» T 4- rn 

rp _ INI 1 IN2 + ±N3 1 IN4 

IN 4 

_ 24.29 + 24.27 + 24.10 + 24.34 

4 

= 24.25 °C 

2. Average Outlet Temperature ( T 0UT ) 

T + T + T + T 

m _ OUT1 OUT 2 OUT3 OUT4 

OUT 4 

39.52 + 39.87 + 37.97 + 41.21 
4 

= 39.64 °C 

3. Temperature Difference in Insulation (AT^g ) 

AT INS = T INS1 " T INS2 “ 49.13 - 38.33 - 10.80 °C 

4. Average Temoerature of Outer Plate (T ) 

T wol + T wo2 _ 64.93+59.19 

wo 2 2 

= 62.06 °C = 335.06 °K 
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5. Average Temperature of Inner Plate (T .) 
— wi 



T . 
wi 



T . 1 + T . - 
wil wi2 



35.34 + 31.66 



= 33.50 °C = 306.50 °K 



Bulk Temperature (T„) 

B 



T + T 
IN OUT 



24.25 + 39.64 



'B 



= 31.95 



7. Mean Temperature Difference (AT) 



AT = T - T = 62.06 - 31.95 

wo B 



30.11 °C 



C. POWER CALCULATIONS 



1. Power Supplied (Q ) 
— P 



V PR V H 
R PR 



(35.580) (3.223) 
2.013 



= 56.97 W 



2. Heat Lost Through Outer Plate (Q ) 



AT INS 

" Lo ax ins /k ins a pl 



10.80 



(0.00635)/(4.18 x 10 2 ) (0.0742)) 



= 5 



.27 W 
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3. Heat Radiated (0„) 

R 



a. Radiation Resistance (R„) 

- .. R 



R 



R 



1 - e 



wo 



1 - £ 



W1 



A £ A__ F . A„ t £ . 

PL wo PL wo-wi PL wi 



i t _i_ + _L. - u = liill 

PL 



A "£ £ 

PL wo wi 



b. Heat Radiated (Q„) 
K 



4 4 

a (T - T . ) 

wo wi 



•R 



R 



R 



(5.669 x 10 8 ) (0. 0742) (335. 1 4 - 306. 5 4 ) 
(1.444) 



= 11.02 W 



4. Heat Convected to the Air (Qair) 
a. Volumetric Flow Rate 



q 



(0.314) (0.4) 
60 



0.00209 m 8 /sec 



b. Density 



p _ 1.0132 x 10 5 

RT (287) (26.97 + 273) 



= 1.177 Kg/m 3 (0.073 lbm/ft 3 ) 
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c. Heat Convected (Q ) 

v ~air 

^air p ^ C P ^ T ouT ~ T IN^ 

= (1.177) (0.00209) (1.0057) (39.64 - 24.25) 

= 0.03825 KJ/sec 

= 38.25 W 

D. AVERAGE HEAT TRANSFER COEFFICIENT 

7- _ ^air _ 38.25 

A ot AT (0.0742) (30.11) 

P-Li 

= 17.12 W/m 2 °C 



E. AVERAGE NUSSELT NUMBER 



Nu 



h d 
K 



air 



(17.12) (0.00635) 
(0.02662) 



4.084 



F . REYNOLDS NUMBER 

1. Mean Velocity (U) 



q _ 0.00209 

A 0.0016 

c 



1.30625 m/sec 
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2 . Reynolds Number ( Re ) 



Re 



pUd = (1. 177) (1.30625) (0.00635) 

y (1. 983 x 10~ 5 ) 



= 492.33 



Re = 492 
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